SPhN-98-01 
CPT-S593-1297 



TIME ORDERING IN 
OFF-DIAGONAL PARTON DISTRIBUTIONS 

Markus Diehl 

CPT, Ecole Polytechnique, 91128 Palaiseau, France 
present address: 
DAPNIA/SPhN, C. E. Saclay, 91191 Gif sur Yvette, France 

and 

Thierry Gousset 

NIKHEF, P. O. Box 41882, 1009 DB Amsterdam, The Netherlands 

present address: 
SUBATECH, B. P. 20722, 44307 Nantes CEDEX 3, France 

Abstract 

We investigate the relevance of time ordering in the definition of off-diagonal 
parton distributions in terms of products of fields. The method we use easily 
allows determination of their support properties and provides a link to their 
interpretation from a parton point of view. It can also readily be applied to 
meson distribution amplitudes. 



1. Recently there has been renewed interest in off-diagonal parton distributions^, i.e. 
in correlation functions of quark or gluon fields between different nucleon states, which 
have been introduced some time ago JIJ. They appear for instance in the description 
of exclusive photon or meson production in 7*p collisions at large photon virtuality Q 2 
and small squared momentum transfer t of the proton [2-5], showing up as the long 
distance nonperturbative quantities when one factorises the transition amplitude into 
short and long distance subprocesses.0 In meson electroproduction there is yet another 
nonperturbative ingredient, namely the meson distribution amplitude which describes 
the transition from the valence quarks to the final vector meson. The corresponding 
diagrams are shown in Fig. |l[ 

The terms nondiagonal, off-forward or nonforward parton distributions are also used in the literature. 
2 The question of factorisation has been discussed in j| for the production of mesons and in Q for 
the production of a real photon. 
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Figure 1: Diagrams for ^*p — > at large Q 2 and small i, where A is (a) a real photon 
or (b) a meson M. They consist of an off-diagonal parton distribution, a hard scattering 
part, and a meson distribution amplitude in case (b). The partons connecting S and H 
have momenta k and k! with plus components fc + = xp + and & /+ = x'p + . 



The relevance of off-diagonal parton distributions is not restricted to these processes. 
In particular the asymmetric gluon distribution has been considered in several 7*p pro- 
cesses at small Bjorken x, such as diffractive production of dijets and of heavy quarks 
|J. For exclusive J/ip production one may relax the requirement of large Q 2 because of 
the large meson mass ]7), as well as for the exclusive production of a Z boson ||. 

2. Let us first have a closer look at the kinematics of the reactions we have in mind. We 
denote particle momenta according to 

l*(q)+p(p) -> A{q')+p(p') , (1) 

where A is a meson, real photon or heavy gauge boson, or a quark- ant iquark pair in a 
diffractive process. One might also have a virtual Z or W instead of the 7* and/or A 
may be charged, replacing the initial or final state proton by a neutron as necessary. We 
will use the momentum transfer A = p — p' from the proton, and the invariants t = A 2 , 
W 2 = (p + q) 2 , Q 2 = -q 2 , m 2 A = q' 2 , m 2 = p 2 . 

We now consider those reference frames where the incoming photon and proton are 
collinear. Picking out any one of them we choose the z axis along the proton momentum 
p and introduce two lightlike vectors 

V = -L(1,0,0,1) , t/ = -±= (1,0, 0,-1) , (2) 

which respectively set a plus and minus direction. We can write 

P " =P V " + 2^ V ' - fr + v " ~ 2 /+(i _ ^ v + A £ ( 3 ) 

and 



« = -=£^, (4) 

having defined the variable ^ = A + /p + , whose exact expression in terms of the invariants 
listed above we will not need here. We further have 

<f = -x NP + v» i + -Q—^ , (5) 
2p + x N 
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where 

x N = . (6) 

l + yjl + ^x 2 B m 2 p /Q 2 

is Nachtmann's and xb = Q 2 /{2p • q) Bjorken's scaling variable. From the positivity of 
particle energies in the final state one obtains 

x N < £ < 1 • (7) 

To ensure that the blob H in Fig. [I] corresponds to a hard scattering we require that 
at least one of Q 2 and m 2 A be large compared with a GeV 2 . For the blob S to be a soft, 
long distance quantity one will need that the transverse bend |A T | is of the order of a 
hadronic scale. 

The threshold region is special in its dynamics: apart from the formation of resonances 
one can expect strong rescattering effects between A and the outgoing proton which 
may destroy factorisation. One will ask for W — m A — m p ^> 1 GeV to exclude this 
region.^ It can be shown that this is equivalent to limiting the range of £ to 1 — £ ^> 
(m p W)/(W 2 + Q 2 ), which according to (f|) also puts an upper limit on —t. Up to 
corrections of order m 2 / (W 2 — m A ) and A^/ iW 2 — m A ) one then has 

Q 2 + m 2 A 

Under these conditions one can find a frame, e.g. the j*p cm., where the initial proton 
is moving fast and collides head on with the photon, p + ^> m p , and where the scattered 
proton is fast as well, p' + ^> m p . Such a frame is natural for a partonic interpretation of 
our process.^ 

3. The nonperturbative transition from the proton to the parton level (see Fig. [1]) is 
described by a two-point function of the form 

M = ^ fdz- e^ +z - (p', a' | T ^(0) 7 + ^) | p,a) , . (9) 

Z7T J z=z V 

This representation holds in the A + = gauge, where is the gluon potential; for 
other gauges one has to insert the standard path ordered exponential between the two 
quark fields. Instead of 7 + in (g) there can be other Dirac matrices, which in order 
to give a leading twist contribution must transform like 7 + under a longitudinal boost. 
Which matrices are relevant depends on the process considered; in virtual Compton 
scattering, 7*p — > jp, one needs for instance 7 + and 7 + 7s- There are also contributions 
from two-point functions where gluon field strengths replace the quark fields in (P). 



3 We remark that we will not actually make use of this condition in the arguments of this paper. 

4 The analysis we perform in this paper may also be done in frames where p and q have nonzero 
transverse components 0, provided that these are sufficiently small, as Ay in the frames considered 
here. 
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We notice that these quark and gluon two-point functions depend on the spins a 
and a' of the initial and final state proton. This spin structure can be made explicit 
by writing (Q) as a sum of Dirac bilinears for the proton multiplied by scalar functions 
Fi(x,£,t). It is these functions that are called off-diagonal parton distributions; we will 
not need their explicit definitions here.Q For ease of writing we will not explicitly display 
the labels a and a' henceforth. 

The usual parton distributions that occur in inclusive deep inelastic scattering involve 
diagonal matrix elements of quark fields or gluon field strengths, i.e. they correspond to 
the limit p' = p and a' = a in @. The variables £ = A + /p + and t = A 2 are then zero 
since A = 0. The diagrams that describe inclusive deep inelastic scattering are obtained 
by cutting the ones of Fig. [I] (a) with the real photon replaced by a virtual one. 

At this point we note that because the proton states in are different the off- 
diagonal distributions do not have any probabilistic interpretation, contrary to the usual 
parton distributions. In particular one cannot expect the Fi to have a definite sign as 
ordinary quark or gluon densities, which are nonnegative for < x < 1. 



4. When calculating the transition amplitude from the diagrams in Fig. [1] the soft blob 
describing the proton-parton coupling translates into the Fourier transformed matrix 
element of a time ordered product of parton operators. We will show that one can 
actually drop the time ordering and instead use ordinary products. 

The same problem has long ago been considered by Jaffe |9|] for the parton distri- 
butions in deep inelastic scattering. For the deep inelastic cross section one needs the 
absorptive part of the forward 7*p — > 7*p amplitude, cutting the diagrams including the 
soft blob, and naturally obtains ordinary products instead of time ordered ones. It was 
however shown in [H that the ordinary product already appears in the full amplitude, 
not only its absorptive part.f] The same result also follows from early work by Landshoff 
and Polkinghorne 0], without being explicitly stated there, and we will apply their 
method to the non-diagonal case here. 

The fact that one can omit time ordering in the parton distributions has important 
consequences for their support properties in the scaling variable x, which in turn are 



crucial for their interpretation in terms of partons || |Tl|]. In the off-diagonal case it is 
also needed if one wants to derive a dispersion relation for the scattering amplitude Q. 
Furthermore it allows one to constrain the two-point function (|9|) using time reversal 
invariance: together with parity invariance one finds that the phases of the scalar parton 
distributions F^x, are a matter of convention, so that they can be defined as real 
valued. 

As a by-product of the demonstration we will obtain the support properties for Fi in 
x and identify different regions in x according to whether the partons are in the initial 



5 A comparison of the various distributions introduced in the recent literature can be found in [Q. 

6 This is also plausible from the point of view of the operator product expansion. There one obtains 
matrix elements of local operators, which are moments of the parton distributions; only the Wilson co- 
efficients describing the hard scattering depend on whether one takes the full amplitude or its imaginary 
part. 
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or final state of the hard scattering HQ This generalises the results obtained in || 
concerning the support and partonic interpretation of the usual densities to the case of 
off-diagonal distributions. 

5. Let us first consider the quark two-point function M of (^|). Our arguments will be 
unchanged if 7 + is replaced by another suitable Dirac matrix. We rewrite 

M = j^yJ d2kTdk ~ Alk+ =*? + (10) 

where 

A= [ d 4 ze ihz (p'\TiP(0)j + iP(z) \p) (11) 



is an amplitude describing the scattering of an off-shell antiquark with momentum —k 
on an on-shell proton with momentum p, 

p(p) + Q(-k) - p(p') + q(-k') . (12) 

Note that this amplitude is not truncated in the parton legs. The plus components 
k + = xp + and k' + = x'p + with x' = x — £ are kept fixed in the following. 

The kinematical invariants A depends on are the virtualities k 2 and k' 2 and the 
Mandelstam variables s = (—k + p) 2 , u = (k + p') 2 and t, the latter being fixed by the 
kinematics of the reaction ([!]). Following |10| we assume that the analytical properties 
of A are the usual ones known from the study of Feynman diagrams, i.e. that it has 
cuts for nonnegative Res, Rev, and singularities for nonnegative Re/c 2 and Kek' 2 . In 
the standard conventions, these singularities are slightly below the real axis of these 
variables.^ 

We can carry out the integration over k~ in ([lO]) taking into account the singularity 
structure of A. To see where the singularities are located in the complex k~ plane we in 
turn express k~ through one of the invariants s, k' 2 , k 2 , u, supplemented by x, kx and 
the variables in (|3|) which are all kept fixed during the k~ integration: 



k' 2 + (k T - A r ) 2 



+ A" 



2p+(x-0 

2p + x ' 
_ u + (k T - A r ) 2 
k ~ 2p+(x+l-Q P 



7 The support of has been obtained in || by analysing the soft blob S in terms of Feynman graphs. 

8 To avoid possible complications due to the dependence of A on the spin of the proton one can 
apply our argument directly to the scalar functions Fi(x, £, t) where the proton spin structure has been 
factored off. The corresponding amplitudes for ( p"2| ) only depend on the above invariants and on the 
plus components p + , A + and k + which are all fixed. The dependence on plus components comes from 
factoring off the proton spin and from the Dirac matrix between the quark fields, and also from the 
choice of gauge A + = 0. 
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There is a sequence of sign changes of the above denominators as x varies from — oo to 
+00. Using < £ < 1 we see that they occur at x = £ — 1, 0, £ and 1. These signs 
determine whether the singularities of A are located above or below the real axis in the 
k~ plane. 

For x < £ — 1 and x > 1 all singularities lie on the same side of the real k~ axis and 
we can close the contour in the other half plane without encircling any singularity. The 
k~ integral vanishes and we deduce that 

M is only nonzero for £ — 1 < x < 1. 

To be able to close the integration contour in the k~ plane we have made the dynamical 
assumption that A vanishes sufficiently fast as k~ tends to infinity in the complex plane 
while k + and are fixed, a limit in which the parton virtualities k 2 and k' 2 become 
infinite. Note that we have not yet integrated over : as remarked in || the integration 
of A over k T and k~ leads to ultraviolet divergences, but this does not contradict our as- 
sumption that the integral over k~ alone is sufficiently well behaved. The full integration 
over k^ and k~ can only be done in a suitably regularised theory; it is known from the 
usual parton distributions that the corresponding renormalisation of the integral (|T0| ) is 
intimately related with their logarithmic evolution in QCD. 

Let us now discuss in turn the regions of x in the interval £ — 1 < x < 1. Details of 
the argument are given in the Appendix. 

A. In the region £ < x < 1, the s singularities are in the upper half plane of k~ whereas 
all other singularities are in the lower half plane. Wrapping the k~ contour around 
the s cuts gives 

r+00 p 

dk'A= dk~ disc s A , (14) 

00 J 

where disc s ^4. is the discontinuity of A across the cuts in s. We use now the standard 
reasoning on the matrix element 



M = — I dz~ e ixp+z (p' I ^(0)7+^(2) I p) 



2tt 
1 



(27TY 



d 2 k T dk~ / d A ze lk - z (p'\ij{0)j + ^{z)\p) 



fc+=xp^ 



(15) 



which does not involve time ordering, inserting a complete set of intermediate states 
between ^(0)7 + and ip(z) as illustrated in Fig. (a). In line with our assumptions 
on its analyticity properties we assume that one can apply the cutting rules to A as 
one would for Feynman diagrams, and obtain 



M = ——^ J d 2 k T dk~ disc s ^l = M . (16) 

In this region of x we can interpret M as the amplitude for the proton emitting a 
quark with light cone fraction x and reabsorbing it with light cone fraction x' . 
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B. In the region £ — 1 < x < 0, the u singularities are alone below the real k axis and 
we can write 

/+oo r 
dk~A = dk~disc u A . (17) 
-oo J 

We now consider the scattering process 

p{p) + q(k') -> + q(k) , (18) 

whose amplitude is given by —A, where the minus sign comes from interchanging the 
quark and antiquark field under the time ordering in (p"T|). Defining 

M' = ~ I dz~ e ixp+z ~ (p' | ip(z) ^(0)7+ \p) . , (19) 



2tt 

and repeating our above argument we find 

M = J d 2 k T dk~ disc u .4 = M' . (20) 

In the present region of x we can interpret M as corresponding to the emission by 
the proton of an antiquark with momentum fraction —x' and its reabsorption with 
momentum fraction —x. 

C. The intermediate region, < x < £, is more involved since the singularities in s and 
k' 2 lie in the upper half plane of k~ while those in u and k 2 lie in the lower half plane. (] 
We notice that in this case x > and x' < 0, i.e. in this situation one has a quark 
and an antiquark flowing out of the initial state. This configuration is specific of 
the off-diagonal matrix element, and in 0, [!| its reminiscence of the quark-antiquark 
distribution amplitude of a meson has been emphasised. 

If we choose to close the k~ contour in the upper half plane we can write the integral 
over k~ of A as a sum of terms due to the singularities in s and in A;' 2 , 

J dk~A = J dk~ ^disc s ^4 + {terms from k' 2 singularities}) , (21) 

where the singularities in k' 2 are cuts and a mass pole, cf. the Appendix. Having 
chosen this representation we consider again M and the scattering (112]), but now 
when we insert a complete set of states, the kinematics allow intermediate states 
that already contain the outgoing proton, \X) = \p'X'). For those states there are 
diagrams which are disconnected at the right-hand side of the cut, cf. Fig. ^| (6), and 
correspond to a cut in the variable k' 2 . They give just the extra terms in ( ^T|) and 
the cutting rules now allow us to write 

1 ( — 
M = J d 2 k T dk~ (disc s .4 + {terms from k' 2 singularities}) = M . (22) 



'The particular singularity structure of A in this region of x has also been remarked in Js], 
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Figure 2: The two alternatives to pick up singularities in the k~ plane, depending on the 
region of x. To the left (right) are the singularities above (below) the real k~ axis. The 
diagrams show the new cut as one changes from one region of x to the next, going from 
top to bottom on the left and from bottom to top on the right. We also give an axis in 
x' to display the symmetry between the two partons. 

In all three regions of x there are of course two ways to pick up singularities in the 
k~ plane: those in the upper half will lead to M and those in the lower half to M'. 
The situation is summarised in Fig. |2|. In regions A and B we have a simple physical 
interpretation if we pick up the s and u cut, respectively, as was already discussed for 
diagonal densities in 0, whereas in region C the choices s, k' 2 and u, k 2 appear as 
symmetric to each other. 

For all values of x we then have 

M M M' , (23) 

i.e. we can drop the time ordering in M and write quark and antiquark field in any order. 
Our argument shows that the origin of this is the integration over k~ , which corresponds 
to fields being evaluated at the same light cone time (z + = 0). We could in fact apply 
it without integrating over at all to obtain the analogue of (|23|) for unintegrated 
parton densities. 
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6. Our preceding discussion can be applied to off-diagonal gluon distributions in an 
analogous way by considering the amplitude for the scattering of an off-shell gluon on the 
proton. Note that one then deals with products of gluon potentials A 11 , while the gluon 
distributions are defined in terms of gluon field strengths . The passage from one to 
the other is however simple in the A + = gauge, since the leading twist distributions 
only involve G +u which in this gauge reduces to d + A v . For a further discussion we refer 
to I. 

In the case of gluons there is no minus sign when one interchanges the order of the 
fields in the amplitude ( |ll|) as there was in point B. above, and correspondingly the 
analogue of M' is defined without a minus sign. 

We also remark that there are symmetry relations for the gluon distributions under 
the change x — > —x' = £ — x in their first argument, so that it is enough to know them 
for x > £/2. In fact the diagrams in Fig. [I] with gluon lines between S and H remain 
the same if one interchanges the light cone fractions x and —x', since H denotes a sum 
of several Feynman graphs. 

7. As mentioned in the beginning, the diagrams of Fig. [1] (b) for exclusive meson pro- 
duction contain not only off-diagonal parton distributions as nonperturbative input, but 
also the meson distribution amplitude, which is expressed in terms of a matrix element 
between the vacuum and the meson state \M). It is easy to show along the lines of 
argument developed so far that the time ordering of quark operators in the distribution 
amplitude can be omitted, just as in parton distributions, as we will briefly outline. 

The distribution amplitude is defined in terms of a two-point function 

— J dz + e - iyq " z+ (M | T 7/5(z) 7 -<0(O) | 0) \ z _ z+ ^ (24) 

or a corresponding one with a different Dirac matrix; in order to suppress the path 
ordered exponential between the quark fields one now needs the A~ = gauge. In 
analogy to ([H]) this two-point function can be rewritten as an integral over It and l + at 
fixed l~ = yq'~ of the amplitude for the transition 

q(l) + q{l') - M(q') (25) 

from the valence quarks to the meson. This amplitude depends on the invariants I 2 and 
I' 2 . Writing 



r 



i 2 + \ 2 T 

2q'-y ' 



/+ = ^r^^- (26) 

and closing the contour of the l + integration in the complex plane one immediately 
obtains that the distribution amplitude vanishes outside the region < y < 1. There 
one can express it in terms of singularities in either I' 2 or I 2 , obtaining matrix elements 
of / 0(z)7 _ '0(O) or — ip(0) %[){z)^~ , respectively. The two possibilities are shown in Fig. |||. 
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Figure 3: The two possibilities to pick up singularities in the / + plane for a meson 
distribution amplitude in the region < y < 1. 

Comparing with the present situation we see again the hybrid nature of off-diagonal 
parton distributions in the region < x < £ : the cuts in s or u are as in ordinary, 
diagonal parton distributions, whereas the singularities in the parton virtualities I 2 or I' 2 
are reminiscent of distribution amplitudes. 

8. The fact that parton distributions and meson distribution amplitudes can be ex- 
pressed in terms of cut amplitudes is important if one wants to show that the s-channel 
discontinuities of the j*p amplitudes shown in Fig. |I] are obtained from appropriate 
cuts of the hard blob H, the soft quantities being "already cut". Already of interest 
in the case of inclusive deep inelastic scattering this is essential for processes involving 
off-diagonal parton distributions, for instance when establishing a dispersion relation [Q. 
We shall not go into details here but wish to make some observations in connection with 
the results we have obtained so far. 

There are different ways to cut the diagrams of Fig. |I| in the overall s-channel, i.e. in 
the variable (p + q) 2 . If the hard blob H is cut in (q + k) 2 , which requires x^ < x < 1 to 
ensure positive energy across the cut, one can cut the soft blob S in s = (p — k) 2 ; if x < £ 
one can also cut S in k' 2 or cut the parton line k' itself in the diagrams. The two-point 
function M can be written in terms of just these cuts. As discussed in the appendix it 
includes in particular the term that comes from cutting the parton line k' and leads to a 
mass pole in the amplitude A of ([y]); hence this mass pole does not appear as an extra 
term in the expression of the discontinuity of the 7*p amplitude. 

H can also be cut in (q — k') 2 , then one must have £ — 1 < x < C, — x^. This comes 
with cuts of S in u = (p + k') 2 , and if x > also with cuts in k 2 as well as the pole 
from cutting the parton line k. In this situation one will use the representation of M 
in terms of singularities in u and k 2 . It is interesting to note that if there is a region 
ijv < x < £ — xn, e.g. in photoproduction of a heavy meson or a Z, one can have cuts of 
H in both (q + k) 2 and (q — k') 2 at the same value of x and needs both representations 
of M at the same time. 

In the case of meson production (Fig. |l] (b)) one can also cut through the blob 
representing the meson distribution function, with one valence quark at either side of 
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the cut, or one can cut through one of the valence quark lines.[^] In this case one will use 
that the meson distribution amplitude can be written in terms of singularities in one of 
the quark virtualities I 2 and I' 2 as discussed in sec. 7. 
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Appendix 

In this appendix we study in some detail the identity of the matrix elements M and M 
introduced respectively in (|l(J) and (|T5|) . As we already mentioned the integration 
has no relevance in the derivation so that we perform our reasoning on the corresponding 
ky unintegrated quantity. We focus on the most complicated region < x < £ where as 
discussed in point C. of sec. 5 the integration over k~ inevitably picks up singularities in 
more than one invariant and where the cutting rules have to be applied with some care. 
In order to deal with complete but readable expressions we ignore here complexities due 
to the spin of both partons and hadrons and introduce 

/+oo r 
dk~ d*ze ikz (p'\T<j)\0)<j)(z)\p} , (27) 
-OO J 

^ r-\-oo r 

M = / dk~ d*ze ik - z {p'\ft{0)<f){z)\p) , 

J — oo J 

where is a charged scalar "quark" field and \p) a scalar "proton" state. 

In the framework we are using we invoke perturbation theory for the purpose of 
analytic properties and of applying the cutting rules. To be consistent we must treat 
quarks as free particles with a mass m q . The scalar analogue of the amplitude A in 
(pT|) then has cuts and a mass pole in each of k 2 and k' 2 . To apply the cutting rules we 
will isolate the poles from the rest of the singularity structure, and also keep track of 
possible disconnected terms of the matrix elements that appear in our argument. To this 
end we use the framework of the reduction formula [HJ (cf. also and introduce an 
interpolating field <P for the proton. For p ^ p' the matrix element (p'\T(fy (0)(fi(z)\p) does 
not have a disconnected part, and in the case of diagonal parton densities one explicitly 

10 Note that such cuts already appear if one takes the hard blob H at Born level, cf. e.g. the diagrams 
in Fig. 2 of @). 
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only takes its connected part. One thus has 



/ ft£ *' ! dWPttMW = ^^-w^-^p, w (28) 

with 

(2n) 4 5 W (p -k-p' + k')-iG=^^- I d 4 zd 4 z' d 4 yd 4 y'e lk - z ~ ik '- z '- tp - y+tp '- y '- 

Z q Z p J 

(□,, + m 2 q ) (D z + m 2 q ) (n yl + m 2 p ) {U y + m\) ■ (0\T<f>\zf)(/>{z) ${y') $\y) |0> , (29) 

where Z q and Z p are the wave function normalisation constants of the quark and the 
proton, and where four- momenta are labelled as in (|12|). In the limit where both parton 
legs go on shell G is a T-matrix element: depending on the signs of k° and k'° the 
partons are in the initial or final state. We emphasise that only the parton mass poles 
in k 2 or k' 2 have been isolated from G which thus still contains the branch cuts in these 
variables, unlike truncated Green functions where the full propagators including these 
cuts are split off. 

To obtain M. we integrate the right hand side of eq. ( p8|) over k~ and close the 
contour in the upper half plane, assuming that convergence is fast enough at infinity. In 
the region < x < £ the only singularities in the upper half k~ plane are due to the 
pole l/{k' 2 — m 2 + ie) and the cuts in s and k' 2 of G, according to our hypothesis on the 
singularity structure. We get 



/+°° / — L \ 
dk~ 27r6(k' 2 - m 2 )JZ q - V iG(pq -> p'q' 
-oo k — Ulq + IE 

/+oo i\t Zq i\ I Zq p „, 

dk~— — — — — — % disc s G + disc^G . (30) 
-oo k — vn £ q — ie k z — + is L J 

Note that the sign of ie in the k' 2 pole has changed from (p8|) to (|30| ) as a consequence 
of separating the contributions of the pole and the cuts, with the pole lying below the 
cuts in the k~ plane. 

Using the cutting rules |T5j (cf. also [III]) the above discontinuities can be expressed 

as 

-idisc s G{pq^p'q') = ^(2tt)V 4) (px + A: / - p') G{pq ^ X)G*{p'q' -> X) ,(31) 

x 

-tdisc k ,iG(pq^p'q') = (2n) 4 5^( Px/ + k')G(pq^p'X')G*(q' ^X') , (32) 

X'^q" 

where px and px> respectively denote the four-momenta of the cut states X and X '. The 
functions G are defined in analogy to (|29|) , i.e. as Fourier transformed Green functions 
with the poles in the external legs removed. They reduce again to T-matrix elements in 
the on-shell limit of the external antiquark legs. The k' 2 cut (|32D has no contribution from 
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a single antiquark state, as expressed in the restriction X' ^ q", where q" has momentum 
—k" : if q' is off-shell such a state is excluded by momentum conservation and if it is 
on-shell then the amplitude G is zero for the non-interacting transition q' — ► q" . 

Let us now turn to the study of hA when < x < £. Inserting a complete set of 
"out" statesQ we obtain 

M = / ^-^(2vr)V 4 )(^ + A;'-p') in (p , |0 t (O)|X) out - out (X|0(O)b) in . (33) 

Due to momentum conservation |X) cannot contain the initial proton state \p) in the 
region < x < £ we are considering. The second matrix element in (|33| ) is then connected 
and can be written in terms of G using again the reduction formula: 

<X|0(O)b> in =— !m_iG(pg-.X) . (34) 
Ar — + 

There are however states that contain the final proton state \p') and lead to disconnected 
parts of (p'\$ (0)\X) . For states \X) = \p"X') with a proton \p") of momentum p" 
we have 

out (/X'|0(O)b') m = 
+ 

If X' is a single antiquark state q", which due to momentum conservation requires q' to 
be on shell, the disconnected term at the r.h.s. of Q3"5D simply reads 

(27r) 3 2p' o 5( 3 )(p'-p") out (g"|0(O)|O) m = (2 7 r) 3 2p' 5( 3 )(p'-p") yfz q . (36) 

For states \X) that do not contain a proton \p") the analogue of ( |34j ) holds. Together 
with the connected term in (]35|) these states give the s-channel cut in A'i, whereas the 
disconnected terms give the cuts and pole in k' 2 . Using 

J2 (p x , + k') = 27r6(k' 2 - m 2 q ) (37) 

X'=q" 

for the pole term we obtain 

11 So far we only had matrix elements between one-particle states or the vacuum, which are equal in 
the "in" and "out" representations, here we need to make the distinction. 



i\ Z„ 



iG(p'q' -> p"X" 



k' 2 



m z + is 



k' 2 — m 2 + is 



iG(q'^X') . (35) 
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M 



{ i — i\ Za 
dk~ { 2n5(k' 2 - m 2 )Jz q y — 



k 2 — m 2 + ie 



iG(pq — > p'q' 



+ E (2n) 4 6^(p' x + k') 

X'^q" 

+ J2W^(Px + k'-p 

X 



i\ Z„ 



—i\Z, 



k 2 — m 2 + ie k' 2 — m 2 — ie 



i\ Z„ 



—i\Z n 



k 2 — m 2 + ie k' 2 — m 2 — ie 



G(pq^X)G*(p'q'^X) 



(38) 



With (0) and the cutting rules (|3T|), (|32|) we finally have M = M, and thus M = M.Q 
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